Loop quantum cosmology and tensor perturbations in the early universe 
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We study the tensor modes of linear metric perturbations within an effective framework of loop 
quantum cosmology. After a review of inverse- volume and holonomy corrections in the background 
equations of motion, we solve the linearized tensor modes equations and extract their spectrum. 
Ignoring holonomy corrections, the tensor spectrum is blue tilted in the near-Planckian superinfla- 
tionary regime and may be observationally disfavoured. However, in this case background dynamics 
is highly nonperturbative, hence the use of standard perturbative techniques may not be very re- 
liable. On the other hand, in the quasi-classical regime the tensor index receives a small negative 
quantum correction, slightly enhancing the standard red tilt in slow-roll inflation. We discuss pos- 
sible interpretations of this correction, which depends on the choice of semiclassical state. 
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I. INTRODUCTION 



One of the major problems of modem theoretical 
physics is how to reconcile the highly successful frame- 
work of quantum field theory with Einstein's general 
relativity, that describes classical gravitational interac- 
tion with astonishing accuracy. There exist several ap- 
proaches in this direction. One, extensively studied by 
the scientific community, is string theory, which attempts 
to unify all forces of nature. Another, dubbed loop quan- 
tum gravity (LQG), focusses on a consistent quantization 
of general relativity alonei'^ii^ii^'^. While the understand- 
ing of physical consequences implied by full LQG is still 
an open research topic, progress has been made by study- 
ing the theory in symmetry-reduced spacetimes. In par- 
ticular, loop quantum cosmology (LQG) opens up a pos- 
sibility of resolving the singularities that plague classical 
cosmological spacetimes^'^'^ i^'^^i^^ . For a spatially flat 
Friedmann-Robertson-Walker (FRW) background with 
a massless scalar field, the model can be analyzed rig- 
orously in terms of physical observablea^^. Although 
the symmetry reduction is performed at classical level 
as in standard mini-superspace quantization, the tech- 
niques used in LQG follow closely those of loop quantum 
gravity. Gonsequently, it can be shown that LQG quan- 
tization is inequivalent to the standard Wheeler-deWitt 
quantization, and thus can lead to significantly different 
predictions^^. 
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The current formulation of LQG mainly deals with the 
quantization of homogeneous spaces but efforts are be- 
ing made to incorporate inhomogeneities at the quantum 
leveli^ as well as within an effective classical framework 
that ensures an anomaly-free constraint algebra even in 
the presence of quantum gravity corrections^. 

In this paper we review the construction of kinematical 
operators and the Hamiltonian constraint on a spatially- 
flat, homogeneous and isotropic background. The ma- 
terial of Sec. |TT] is a self-contained review of the litera- 
ture, but also implements, with the inhomogeneous case 
in mind, the p. ^ p^" quantizatior^i^ both in the grav- 
itational and matter sectors. In Sec. IIIII we use recent 
flndings on cosmological perturbations to compute phys- 
ically relevant quantities of the early universe. In partic- 
ular, we concentrate on tensor perturbations originated 
at an early stage in the cosmic evolution in the pres- 
ence of inverse- volume corrections (holonomy corrections 
ignored). The near-Planckian regime of LQG typically 
leads to a strongly blue-tilted tensor spectrum, which is 
disfavoured by observations if the tensor-to-scalar ratio 
is not too small. However, in this regime the dynamics is 
intrinsically nonperturbative, as it takes place near the 
bounce, and the use of cosmological perturbation theory 
may not be completely reliable. There are also a number 
of other assumptions, made clear during the discussion, 
which may be relaxed in a more complete (and compli- 
cated) analysis. On the other hand, the quasi-classical 
regime of LQG is close enough to general-relativity re- 
sults to be in agreement with data. In particular, the ten- 
sor index receives a small negative quantum correction, 
slightly enhancing the standard red tilt in slow-roll infla- 
tion. Both regimes are concrete realizations of effective 
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dynamics with modified dispersion relations, which have 
been long argued to play a role in the trans-Planckian 
problem for inflatio n^^'^^i^^ . 



II. INVERSE- VOLUME AND HOLONOMY 
CORRECTIONS 

A. Classical canonical variables, volume, and 
gravitational Hamiltonian 

In Ashtekar-Barbero formulation of gravity, the canon- 
ically conjugate variables are the connection and the 
densitizcd triad liS, where a,i = 1,2,3. Here a is 
the spatial index whereas i denotes the directions in the 
internal space. Being conjugate variables, and Ef 
satisfy the following Poisson bracket 



(1) 



where G is Newton's constant of gravitation. The physi- 
cal volume can be expressed as 



V = y" (fxVdetE, 



(2) 



where det E 



eabce''''EfE]El/i\. eabc is the Levi- 
Civita symbol: €123 = £231 = £312 = +I7 antisymmet- 
ric in pairs of indices. Repeated upper-lower indices are 
conventionally summed over. 

The spatially-flat Friedmann-Robertson- Walker met- 
ric is ds^ = ~dt^ + a{t)^dx^. A physical distance d is 
determined by the dimensionless scale factor a and the 
coordinate (comoving) length xq, d = axo- The volume 
of a spatially-flat universe is inflnite and thus, to render 
the volume integral ^ well deflned, we need to consider 
a patch of the universe with finite 'fiducial' coordinate 
volume, say Vo- The physical volume of the patch then 
reads V = a^Vg^. 

In this background, the symmetric Ashtekar connec- 
tion is expressed as A^^ = cVq^^^ujI^, where cVq^^^ = 7a 
classically, 7 is the Barbero-Immirzi parameter— i^, and 
a dot is derivative with respect to synchronous time 
t. Here are the components of invariant 1-forms 
that satisfy w^cj^Jy = Sab- The extrinsic curvature 
is i^^ = (A^ — rj,)/7. Here T^ is the spin connec- 
tion, which vanishes for the spatially flat FRW space- 
time. The symmetry-reduced conjugate momentum is 



E? 



pVq ^^'^e°, where pVq 



(the sign of p is cho- 



sen to be positive by convention) and ef (often denoted 
°ef) is the invariant vector field dual to w^: efcol = S^. 

In the quantum theory the connection itself cannot be 
promoted as a well-defined operator; only the holonomy 
of the connection can be consistently quantized. The 
holonomy along an oriented path e is an element of the 
SU{2) group defined as 



he = V exp 



dXe%X)A''^iX)Tk 



(3) 



where dx"" /d\ = e°(A) is the tangent vector along the 
path parametrized by A and Ti, i = 1,2,3, are the 
three generators of the su(2) algebra in irreducible j- 
representation (of dimension 2j + 1). 

For technical convenience, in the quantization of the 
scalar constraint operator we fix the representation to 
be the fundamental one, i.e., j = 1/2 (this choice may 
be justified also by theoretical considerations^!'-^, which 
however are not compelling and can be bypassed). In 
this representation, the generators are — crfc/(2i), <7k 
being Pauli matrices, and satisfy 



TiTk 



4 ' 



(4) 



where I is the 2x2 identity matrix. In the FRW case, 
the coordinate space is flat and the tangent vector is the 
same at every point along the edge e. In particular, the 
holonomy of oriented length £q along the i-th direction is 



exp(£oVo ^^^cn) = exp(/icTj) . 



(5) 



where the path ordering V becomes trivial and p = 

1/3 

Iq/Vq is the ratio between the holonomy length and 
the size of the coordinate volume. To keep notation light, 
from now on we omit the superscript (/i) in the holonomy 
symbol. 

In view of implementing lattice refining quantization, 
it is convenient to define a new pair of variables a a^^'^'^ 



2p 



(6) 



where p, = fj,{p) is an arbitrary dimensionless function of 
p. Later arguments will suggest 



'A 

= I - 

P 



PI 



2n 



0<7i<l/2, (7) 



where A is a constant of dimension (length)^ which can 
be coordinate dependent, and = Gh is the squared 
Planck length. Then it is easy to show that b and w are 
also canonically conjugate and 



{5, ^}^!^(i + „) 



(8) 



Defining A = /i//i as the path length in units of /i, we can 
express the holonomy ^ in terms of the new variables 
as 



hi — exp {2XbTi ) = I cos Xb + 2Ti sin Xb , 



(9) 



where we employed the identity Eq. (jj]). In quantizing 
inverse powers of the scale factor, one faces the problem 
that naive inverse operators (p)"", s > 0, fail to remain 
well-defined. As a standard procedure, one rewrites in- 
verse powers of the densitized triad using classical Pois- 
son brackets involving positive powers of p. In particular, 
the following classical relation is useful: 



{6, w'} 



iirGj 



1 -I- n)lw 



i-i 



(10) 
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valid for arbitrary values of /. Via Eqs. ([9|) and (jlOp. one 
can express powers of w as 



..i-i 



j:MnK{h-\w'}) 



SirG-fil + n)lX 
3 



cos Xb {sin Xb, w'} 



87rG7/(l + n)A 

— sin Xb {cos Xb, w'}] , 



(11) 



(12) 



where we wrote explicit sums when summation conven- 
tion is not clear. To get this result, we used the algebraic 
relation ^^TiTi = — C2(j)I and took the trace in the 
j = 1/2 representation (for which the quadratic Casimir 
invariant C2 ~ 3/4). 

For the purpose of quantization, the Hamiltonian con- 
straint needs to be expressed in terms of holonomies. We 
review this calculation^ i^^'^^i^^ , starting from the expres- 
sion of the classical Einstein-Hilbert Hamiltonian con- 
straint IGTTGHg = J d^KNCg-. 



Ca 



VdetE 
1 



7^ VdetE 



ah 



6 



7'Vo 



(13) 
(14) 



where N is the lapse function (the choice N = \ corre- 
sponds to synchronous time t, whereas = a to con- 
formal time r), Fl^ = daAl - dbA^ + e)^AiA^ is the 
gravitational field strength, and indices in square brack- 
ets are antisymmetrized. In the second and third equali- 
ties we have made use of maximal symmetry of flat FRW 
backgrounds. 

To write the field strength F^j, in terms of holonomies, 
one considers the holonomy /iq^j. = hjhkhj^h^^ on a 
closed oriented square path whose edges have length /i 
and are labelled by spin indices j, k, j, k. Ther^i^i^li^^ 



tr(T,/in^.J 



(15) 



and TiTjTi = 
'^Sijk sin^ 2Xb. 



Tj/4, one can show that 



Using Eq. ^ 
tr(Ti/in^., ) 

The terms in Eq. (|13p involving the densitized triad 
E°' can be manipulated as follows. Starting from Eq. ^ 
and the property eibcf-"^^^ — 2(5j?jj(5^j, one can prove the 
following formula due to Thiemann: 



rpa rpb 



Vd^ SEf 



2ttGj£o 

(16) 

where in the second equality we used the expression 
hkSh'i^^ /SA"^ = —ioe'j.Ti. In terms of the canonical pair 



(6, w), 



7fc ^J^fc 



^abc 



VdctE 



/^n \ 2(l + „ 



2 J 2TTG-f£o 
cos Xb {sin Xb, w^ci+'i) } 

- sinAfe{cosAfe, w^nr^} ^(17) 



in agreement with Eq. (fT2| . The total Hamiltonian con- 
straint reads 



^ i.k' 

'0 ^ 



= lim 



^0^0 nG-f^el \ 2 



— sin^ 2A& 



sin A&{cos A&, } 
- cos A6 {sin Xb, ^ W^^T } 



B. Quantization and lattice refinement 



(18) 



In the quantization procedure of loop quantum cos- 
mology, the classical Poisson bracket is replaced by the 
commutator bracket as in standard quantum mechanics, 



(19) 



As mentioned earlier, the connection-dependent operator 
6 does not have a well-defined action on the kinematical 
Hilbert space. Nonetheless, the other elementary variable 
w can be promoted to a self-adjoint operator. 



w^w = — - — (1 + n)— . 

61 db 



(20) 



It is easy to check that the functional [b) = e*'*''' = 
(6 1 A) are the eigenfunctions of w. In Dirac notation, we 
can write down the action of Eq. (1^ on the states |A) 



as 



w\X) 



(1 +?i)A|A) = w\X) 



(21) 



These eigenstates form a basis {|A)} of the kinematical 
Hilbert space. On this basis, the action of a holonomy 
operator of edge A' can be constructed from the basic 
operator e*^''', which acts simply as a translation: 



s'^'^IA) = |A + A'). 



(22) 



We are interested in computing the spectrum of well- 
defined inverse- volume operators, so we consider the val- 
ues of I being < Z < 1 in Eq. (fT^ . The ambiguity 
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parameter / determines the initial slope of the effective 
geometrical density. To preserve coordinate invariance 
when quantizing geometrical densities before symmetry 
reduction, / must be discrete, Ik = 1- (2fc)~\ k e N^^'^". 
Hence one can select the bound 1/2 < ^ < 1, which is 
also favoured phenomcnologically^S; a natural choice is 
I = 3/422^. 

To quantize Eq. (|18p in terms of elementary operators, 
we fix the length of the holonomy to be unity (A = 1), 
and consider their symmetric ordering so that Eq. 
becomes the self-adjoint operator 



li-i _ 



3i 



87r4j7/(l + n) 



cosfe |w|'sin6 — sin 6 |w|'cos6 



.(23) 

The absolute value of w is taken in order for the eigen- 
values of to be real. In fact, it is easy to check 
that the basis states |A) are also the eigenstates of the 
operator Eq. 



87r4i7(l-Hn) 



x(|A+l|'-|A-l|0|A). (24) 

Although the discussion so far has been at the kinemat- 
ical level, the motivations for taking A = 1 and ft cx 
are mainly dynamical. They will be the subject of the 
remainder of this Subsection. 

If translation invariance is broken, e.g. by inhomoge- 
neous perturbations, the comoving volume of the system 
under inspection can be discretized as a lattice whose J\f 
cells have characteristic comoving size io and correspond 
to the vertices of the spin network associated to Vo- We 
identify the ratio of the cell-to-lattice size io = {Vo/f^Y^^ 
with the previously ad-hoc function p.(j)), under the re- 
quirement that the lattice be refined in time. Hence 

The rescaling of the fiducial volume sometimes led to 
the concern that LQC quantum corrections (as well as all 
statements regarding the scale at which transition from 
semiclassical to quantum regimes occur) break confor- 
mal invariance of a fiat FRW background; scaling invari- 
ance can be realized if proper Vo factors are considered 
in the observables but then the theory and its observ- 
ables would depend on the choice of the fiducial vol- 
ume. However, in the presence of inhomogeneities there 
is no conformal freedom and, on the other hand, fiuxes 
through two-surfaces F — cPaEln'^Ti [ua being the 
normal to a surface S), that is, elementary areas, are de- 
termined by the inhomogeneous spin-network quantum 
state of the full theory associated to a given patch"^^ . This 
implies that to change the fiducial volume a^Va would 
change the number of vertices of the underlying physi- 
cal state. Therefore there is no scaling ambiguity in the 
equations of motion^i^^, although the physical observ- 
ables (through A, see below) will depend on the choice of 
spin-network state. 

A feature of loop quantum gravity, which is not an 
assumption but a consequence of the full theory, is that 



the spectrum of the area operator is bounded from below 
by the Planck scale. Although the area spectrum of LQC 
does not share this property, one may draw inspiration 
from the full theory and consider lattice cells to be such 
that their kinematical area is bounded as 



[alaf > Api = 2V37r7# 



i2 

PI 1 



(25) 



where the numerical constant stems from the direct cal- 
culation in LQG'^-'^'^. Big-bang nucleosynthesis can place 
a bound over the smallest physical shown in 

Ref.'^^. If the inequality is saturated (smallest possible 
holonomy path), the comoving cell area 



oc 



1/3 



(26) 



is also the comoving area gap, that is, the smallest nonva- 
nishing eigenvalue of the area operator measuring comov- 
ing surface a^^'^'^1'^^ . As the Universe expands the comov- 
ing area gap shrinks to zero and the geometry is better 
and better described by classical general relativity, while 
near the big bang quantum effects become important. 
However, there may be more general forms than Eq. (j26p . 
as the spin-network state described by the lattice can be 
(and usually is) excited by the action of the Hamiltonian 
operator on the spin vertices, increasing their number 
and changing their edge label o^^'^^ . This process has not 
yet been established univocally in the full theory, so it is 
convenient to parametrize the number of vertices as^^ 



U{t) = u{t)v/il„ 



(27) 



where the (dimensionless) function f^{t) is state de- 
pendent (the density of vertices per Planck volume) 
and, by assumption, coordinate independent; its time- 
dependence is inherited from the state itself. As the 
kinematical Hilbert space is usually factorized into gravi- 
tational and matter sectors, the problem here emerges of 
how to define a natural clock when matter does not enter 
in the definition of a (purely geometrical) spin network. 
This issue will require a much deeper understanding of 
the theory; so, as unsatisfactory as Eq. (P7|) may be, we 
take it as a phenomenological ingredient in the present 
formulation of inhomogeneous LQC. 

From the definition of ft, one gets Eq. ([7]) if ftp{t) scales 
as a^^^"~^)(t). A is some physical squared length deter- 
mined by the theory. Even in the improved quantization 
scheme, it may differ from the mass gap Api by some 
0(1) constant and, as there is not yet a general consensus 
on this issue, we shall set A = Api only in Eqs. ([571) - ([M)) . 

Homogeneous models adopting Eq. (|26|) feature 
holonomies which depend on triad variables; in other 
words, curvature components are constrained by the area 
operator although this does not appear in the full con- 
straint. On the other hand, in inhomogeneous models 
the dependence of the parameter ft on p is implemented 
at state (rather than operatorial) level, in closer confor- 
mity with the full theory^'^. For these reasons, we will 
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keep n general within a reasonable interval which can be 
argue to be < n < 1/2. In fact, if = 1 (n = 0) the 
states I A) = coincide with the basis eigenstates of the 
momentum operator p, with eigenvalues w oc piSii^. The 
lower bound n — Q comes from the definition of p, (the 
number of vertices must grow with the physical volume) 
and corresponds to a configuration where the total num- 
ber of vertices is kept constant during the evolution of 
the patch. This case leads to severe restrictions of the 
matter sector if the wavefunctions solving the Hamilto- 
nian constraint are required to be normalizable and to 
reproduce the classical limit at large scales^; so it will 
be excluded. In the 'improved quantization scheme — 
n = 1/2, {|A)} become the eigenstate basis of the volume 
operator, w oc p^^"^ — V. In this case, the number of 
vertices grows linearly with the volume (constant fluxes, 
constant spin labels along lattice edges). There can be 
configurations with n > 1/2 but they would correspond 
to states with high and decreasing spin labels, which are 
not realized generically. 

To calculate the curvature at the lattice sites within 
Vo, we need to specify closed holonomy paths around 
such points. A generic holonomy is given by the composi- 
tion of elementary holonomies over individual plaquettes. 
Therefore we set the length of the elementary holonomy 
to be that of the characteristic lattice cell, n — fl. 

In this case, the quantum Hamiltonian (self-adjoint) 
operator corresponding to Eq. is wcU-dcfincd: 



Finally, we define the function 
27^ A' 



a 




l+4n 

2(1 + ") 



A, 



(33) 



which is nothing but a correction function to the eigen- 
value of the operator p2/p3/2 (quantization of Eq. ([T7)) ') 

on |A), p2/p3/2|A) = ay/p\X). Note that, a priori, another 
legitimate choice is to write Eq. (|T6)) as 



Vpoc — cx —— 

r op riQ 



^tr(r,/i,{/i-i, V}), (34) 



which introduces another ambiguity < r < 1. Then it 
is easy to show that 



1 -I- n , 
3r 



1 



1 



(35) 



C. Effective background dynamics 



VoCg = -4 sin 25 A sin 26, 



(28) 



where 



A = 



3i 



A" 



327r^2^73 V 2 

X cos6 |w| 2(1+") sin& 



w 1+" 



sin6|wh(i+")cos6 



and we used the relations 



A" \ ~ 

w 

2 



/2 = 2 



A" \ ~ 



2w" 



(29) 



(30) 



Acting on a state |A), which is an eigenstate of A with 
eigenvalue Ax, one has 

VoCg\X) ^ Ax+2\X + 4:)-iAx+2+Ax-2)\X)+Ax-2\X-4:) , 

(31) 



where 



A, 



1 + n / A" \ 2(i+„) 
"8^ 



8^4i7(l + n) 



2(l + n) 



3(1 + 2,1) 

xA 2(1+") 



1 - 



(32) 



Effective cosmological equations of motion are derived 
from the expression of the effective Hamiltonian con- 
straint. To get the latter, typically one considers a semi- 
classical state IV'sc) peaked around some point {w, b) in 
the classical phase space. In this state one then computes 
the expectation value of the Hamiltonian constraint op- 
erator using appropriate inner product. Accordingly, for 
the gravitational Hamiltonian operator Hg we may write 



(V'scl-ffgl^As. 



3A^ ^ sin 



-2 9 



(36) 



where we have approximated the expectation value of 

the operators in Eq. (j28p as (■(/'sc|sin/2c Asin/ic|V'sc) ~ 
Ax sin^ p,c. 

Including matter components, imposition of the effec- 
tive Hamiltonian constraint (that is, variation of the total 
effective Hamiltonian i?tot — {Hg) + (ffmat) with respect 
to N) yields 



■ 2 - 

sm /ic 



SttG 



■P, 



(37) 



where p — p ' 5{H-aiax) /5N . When a — 1 and the mat- 
ter sector contains only a massless free scalar field (po- 
tential V{(p) — 0), then Eq. ([37)1 is exact^. In general, 
however, the evolution of a finitely-spread semiclassical 
state will produce quantum fluctuations leading to ad- 
ditional corrections to Eq. (|37p . Here, for simplicity, we 
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shall ignore such contributions; these were discussed in 
Refs42.iii. The Hamilton equation of motion for the den- 
sitized triad p = {p, -fftot} gives the Hubble parameter 



H 



a sin 2/zc 
— = a — — — 



(38) 



In the classical limit, c 'ya. Combining Eqs. ((37)) and 
one gets the Friedmann equation 



8ttG 



P_ 

Pc 



where 



SttGj^^'^P 

From the definition of the number of cells/vertices, 



2/3 

Gpc oc ( — 



(39) 



(40) 



(41) 



The improved quantization scheme (— , Eq. (^5)1 ) is the 
only one where pc is constan t-'^^i^^'^'^ , as A/" oc V. Other- 
wise, it will depend on the function /^(i) wherein we have 
encoded our ignorance about the underlying state''^. In 
both cases, the critical density is a number density which 
does not depend on the size of the fiducial volume, nor 
on coordinates. 

We assume validity of the semiclassical approximation 
and that the semiclassical wave packet of the universe 
does not spread appreciably, so we presently stick to 
Eq. also in the presence of a nontrivial potential. 
For a scalar field Lp with momentum -k^ and potential V , 
matter Hamiltonian is 



N 



P 



which upon quantization can be written as 
'1 



N ( -p-3/27r2 + p3/2V 



(42) 



(43) 



The effective matter Hamiltonian can be viewed as the 
expectation value of the matter Hamiltonian operator in 
the semiclassical state {tpsc) (whereon acts multiplica- 
tively): 



(V'scl-ffmatlV'sc) « N 



2p3/2 



-P' 



3/2y 



(44) 



where we have neglected quantum fluctuations from 
the matter sector and h' = (A|p3/2p-3/2|^^ _ 
((A|^Ji^i^|A))3/[2(i-0(i+")]. Using Eq. 



(45) 







I 






\-( 
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21 y 











2(l-i)(l+n) 



The equation of motion ip = {(p,Htot} yields tt^ 
p^/'^ip/{Niy), while ir^ = {7r<p,i?tot} leads to 



0. 



(46) 



As > has a maximum at A = 1 and then decreases 
down to unity for large A, the friction term in Eq. (j46p 
changes sign during the evolution of the universe, the 
flrst stage being of superacceleration, as anticipated. 

Since the momentum operator in w — 2p//I, the total 
p-dependencc of A is 



A 



3V3 



P\ 



l+n 



2(1 + n) Api VA/ 



(47) 



Then the eigenvalues of the volume operator V are pro- 
portional to A3/(2+2n)^ and the classical limit (large- 
volume approximation) corresponds to A 3> 1. Consis- 
tently, in this regime the eigenvalues of jwl'^^ are « n/"^ . 
Then Eq. psp can be approximated as 




3r 



2(1 + n) 



1 
6A2 



l + a 



where 



and 



ac 



4(1 + n). 



[3r - 4(1 + n)] [3r - 2(1 -I- n)] f Api 

A~ 



3^2 

Near Planck scale (A <C 1 



a w A^ 



where 
g„ = 4(1 -f n) - 3r , 

On the other hand, when A ^ 1, 
2-1 



aV. 



1/3 \ 9° 



(48) 



(49) 



(50) 



(51) 



3\/3 A 
2(1-Hn) A^ 



2(l-hn)A^ 



= l + Vc 



aV, 



1/3 



where 



2(l + n)(2-0 / Api 
27 V A 



Near Planck scale (A ^ 1) 



3(2-0 

V K. A2(l-i)(l + ") = J/„ 



aV, 



1/3 \ 9- 



(52) 
(53) 

(54) 
(55) 
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where 



i{2-l) 
l-l ' 



3V3 A ' 
2 + 2n Api 



2(1 + ") 



Assuming A = Api, the parameter ranges are 



4 < c < 6, 

-0.01 0.11, 

33/4 27 
1.61 w — ^ 6.75, 

\/2 ^ 4 



1 < < 6 , 
2 

0.07 « — < 
27 

Vq > 3^/2 « 5.20, q^>9. 



0.07 « J- < < ^ « 0.17 
27 6 



(56) 



(57) 
(58) 

(59) 
(60) 
(61) 
(62) 



Although one can resort to different quantization 
schemes, general arguments show that Eqs. pS)) . ([5T|). 
((53)) . and (|55p maintain the same structure, where the 
coefficients c, Qq,, and are robust in the choice of the 
parameters, inasmuch as their order of magnitude does 
not change appreciably 4-. AH these parameters can be 
set to their 'natural' values, which are dictated by the 
form of the Hamiltonian or other considerations. From 
the calculation leading to a and j/, it is clear that the 
natural choice is 



/ = 3/4, 
which gives c = 6 and 

36 



1, 



1/2, 



15.59, 



(63) 

= 3 , 
gi/ 15 . 

(64) 



At this point it may be useful to summarize different 
parameter choices in the literature: the same sequence 
of steps we reviewed was followed in Refs.— (arbitrary 
j, I = 1/2, n = 0; a = l)r^ (arbitrary j and I, n = 0; 
a = 1)P {j = 1/2,1^ 1/2, r^l,n= 1/2)^ (arbitrary 
j, I, r — 1, n — 1/2), and partly'^* {j — 1/2, arbitrary I, 
r = 1 , arbitrary n) . 

All quantum corrections depend on the ratio 



£ 



PI 



Vl/3 



1/3 



(65) 



In a purely homogeneous universe, the left-hand side is 
volume dependent and inverse-volume corrections will 
feature V, even in the case n = 1/2 (when r = 1, a = 1 
but in general r ^ 1; v is never constant for any al- 
lowed value of the parameters). For a compact universe 
(e.g., a sphere or a torus) the fiducial volume is iden- 
tified with the total physical volume, which is given by 
the theory; so no issue arises in this case. Otherwise, 



the fiducial volume is arbitrary and hence unphysical. 
However, in the lattice-refinement picture the same quan- 
tity is determined by the inhomogeneous state through 
the number of vertices Af and the function /^,. Inho- 
mogeneity is not simply invoked to justify Eq. (|65p but 
can be implemented via a spatially local metric as in the 
covarianf^-^iiSiiiSi^i^ and separate universe approach^^ to 
cosmological nonlinear perturbations. 



D. Background expansion 

In order to solve analytically the perturbed equations 
of motion of the next Section, it will be necessary to 
choose a simple background. In standard inflation, this 
is achieved either in de Sitter spacetime (a = e^*, H 
constant) or when the scale factor has a power-law be- 
haviour: 



a{t) (xtP, p>0. 



(66) 



We should mention here that a is singular as t ^ 0; thus, 
close to bounce point Eq. fails. We comment on this 
issue in the discussion. 

In conformal time t = J dt/a Eq. (|66p would be a cx 
-rP/(i-p) ifp^ 1, but for simplicity we define p — p/ (1— p) 
and make the ansatz (up to an arbitrary normalization 
costant) 

a = TP. (67) 

Below, the symbol p should not to be confused with the 
triad, which will be always written as a^. Then H = 
drCi/a = aH ~ p/t. The first slow-roll parameter reads 



H 



= 1 



1 
P 



(68) 



Inflation occurs for p < —1, superinflation when — 1 < 
p < 0, and normal expansion when p > 0. In particular. 



1 



P 



(69) 



In conformal time, the de Sitter solution {p — —1) is 
a = {—Ht)^^, H = — l/r. de Sitter and power-law so- 
lutions can be used as sensible backgrounds only if they 
are classically stable against homogeneous perturbations. 
This check was done in Ref.— . 



III. TENSOR PERTURBATIONS AND 
OBSERVABLES 

Metric perturbations can be decomposed into three dif- 
ferent modes: scalar, vector and tensor. At linearized 
level, these modes do not couple to each other and gauge 
transformations do not mix them among themselves. 
Thus, they can be studied independently. 
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In considering quantum corrections to the effective 
Hamiltonian constraint, one needs to ensure that the 
modified constraint algebra remains anomaly-free in the 
canonical approach. In classical theory, closure of the 
constraint algebra is guaranteed by general covariance. 
However, in the presence of quantum corrections, one 
needs to explicitly show that the constraint algebra can 
be made anomaly-free. Such analyses have been recently 
performed in the presence of inhomogeneous perturba- 
tions including only inverse-volume corrections^^. 

Due to inverse-volume and holonomy corrections, it 
was shown that vector perturbations are suppressed even 
faster than in classical cosmology^. The derivation of 
gauge-invariant variables for the scalar sector are still 
under investigation. So here we focus only on tensor- 
mode dynamics with inverse-volume corrections. From 
now on we fix Vo = 1; volume prefactors will be restored 
at the end of the calculation. 

The linearized equation of motion for tensor modes 
has been computed^. There, triad and connection com- 
ponents are separated into a FRW background and an 
inhomogeneous perturbation. 



2 ra 



One can see that 



5Ef 



SKI 



c5l + ^Tl+^5K). (70) 



-drK 



-hf 



(71) 



where ft," is the transverse traceless part of the perturbed 
3-metric and the curvature and triad perturbations are 
canonically conjugate: 

{SK{^),SE'^{y)} = 87rGSis;S{^,y) . (72) 

One can go to momentum space and denote with hk ei- 
ther of the two polarization modes of the graviton h^, 
where k is the comoving wave number of the perturba- 
tion. When only inverse-volume corrections are taken 
into account and in the absence of anisotropic stress, the 
equation of motion for hk is^ 



d^hk + n 2 



din a 
din a 



drhk + a^k^hk=0. (73) 



Near Planck scale 



2gc 



9^ftfc+H(2-<7,)a,ftfe+a^^ k'hk^O, (74) 



drhk 



while in the quasi-classical regime 



dlhk + n 



2 + ca 




l + 2ae — k'hk^O, (75) 



where in the last equation we have retained only leading- 
order terms in a^. 



A. Near-Planck regime 



To solve Eq. ((74)) . we make some field and variable 
redefinitions. First, we define the Mukhanov field Wk = 
ahk, the constant 1/ = l/2—p/(l+pqc,) (where p^Q, 7^ — 1), 
and recast Eq. ([74)) in the variable z = J dra = ra/(l -I- 
PQa)- Then 



4i/^ 



1 



4z^ 



Wfc = 0, 



(76) 



and the solution reads (Ref.— formula 8.491.5) 

Wk = CiV^.Hi'^\-kz) + C2V^^Hj^H-kz), (77) 

where Hi, are Hankel functions of the first and second 
kind. Some useful formulae to get Eq. ([76]) are: 



drhk = - {drWk ~ Ti-Wk) , 

a 

dlhk = - [dlwk - 2ndrWk + eH^Wfe] 



dr = ad^ , 



(78) 

(79) 
(80) 



In the classical limit, z ^ t and one recovers the usual 
result. Taking C2 = (in order to get the correct asymp- 
totic behaviour at small scales, advancing plane wave), 
one can study the large- and short-wavelength limits of 
the solution {i' > 0) 



Wk 



Wk 



\kz\ < 1 , 



Ci^e-<'''-+^''+i) , |fcz|>l. 



The tensor spectrum is conventionally defined as 
Vh k^ 



100 2007r2a2 



+ .X 



k=k. 



(81) 
(82) 

(83) 



where k.^, is the wavenumber at horizon crossing, the solu- 
tion is at long wavelengths, the sum is over the two polar- 
ization modes, and angular brackets denote the vacuum 
expectation value of the operator Uk = ahk- The con- 
stant Ci is determined by choosing the Bunch-Davis vac- 
uum (asymptotically Minkowski) in the classical regime, 

Wk ~ e-^'^Vv^, as follows. Taking Eqs. jTll), (Ell), and 
()19p . the commutation relation between Uk and its con- 
formal derivative reads 



[Ufci, drUk^] = 327r^|iiQ!(5(fcl, fc2) . 



(84) 



The Mukhanov variable is expanded in terms of creation 
and annihilation operators, 



ilk = Wkttk + w^al, 



(85) 



where * indicates complex conjugate, Wk is the classical 
solution, and 



4J = S{ki,k2), 



(86) 
(87) 
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Plugging equation (|85p into ([84]) and using equation ([86]) . 
one gets 



Hence, from Eqs. ([Ml) and (ESI) |Ci| = y/8TT^£'f,i/k. 

Horizon crossing is defined at the moment when per- 
turbations are frozen, namely. 



1 



2z 



n 

a 



1 

P 



(89) 



In the classical limit, fc* ^/2aH, as expected. The 
above expression is well defined only ifp> 1/(1 — g^), 
which, according to the value of q^, will determine (com- 
patibly with the stability of the background solution) 
whether the universe superaccelerates or not. In order 
to avoid interpretational issues^^, the stronger condition 
p > —1/qa must hold (superinflation, if the universe ac- 
celerates), so that the time variable z flows along the 
same direction at r and modes exit (rather than enter) 
the Hubble horizon at crossing. Then, 



^ l^2{l+p+pq^)/(l+pq^) ^ 



and the tensor spectral index is 



din 



din fc 



k=k. 



2(e + ga) 
e + qa-1' 



(90) 



(91) 



where / = fc/(27r) is the frequency of the signal, /o = 
aoHo/{2n) ^ 3.10 x lO'^^ Hz, h = 0.716, and'r = 
A^/Ag is the ratio between tensor and scalar (^1) am- 
plitudes. Observations of pulsar timing, interferometer 
experiments (LIGO, LISA), and the theory of big-bang 
nucleosynthesis (BBN) can place strong constraints on 
the tensor index^^. For instance, taking the upper bound 
r < 0.30, from pulsar timing nx ^ 0.79, while from BBN 
riT < 0.15. These values are incompatible with the above 
predictions of near-Planckian LQC. Therefore, modulo 
some important issues we shall comment later, we might 
conclude that a near-Planckian accelerating phase might 
have occurred only at very early times (unobservably 
large scales), and for a short period. 

In homogeneous models the duration of this regime 
was held to depend on the spin representation of the 
holonomies, small j implying a very short superinflation- 
ary period^^. In inhomogeneous situations, this problem 
is reinterpreted and relaxed in terms of the lattice em- 
bedding of Ref ji^. The volume spectrum depends on the 
quadratic Casimir in j representation: ~ y'^/^ ~ 

y/Ozij) ^ j- A higher- j effect can be obtained as a 
refinement^'^ of the lattice (smaller /x), thus allowing for 
long enough superacceleration. A change in p,{p) can be 
achieved by varying the comoving volume Vq. This is 
an arbitrary operation in pure FRW, while in inhomoge- 
neous models /I is a physical quantity related to the num- 
ber of vertices of the underlying reduced spin-network 
state. As long as a calculation of this effect from the 
full theory is lacking, we will not be able to predict the 
duration of the small- volume regime. 



In the quasi-de Sitter limit p « — 1, ht ~ 2qa/{qo, ~ 1), 
leading to a blue-tilted spectrum riT > 12/5. The tensor 
index is greater than 1 also when the universe is in deep 
superacceleration (e ^ —qa, P 0~), ut ~ 2. Prior to 
the calculation of tensor perturbations, some toy models 
already favoured a strong blue tilt of the spectrum^i^. 
This state of affairs is typical of LQC inflation and su- 
perinflation near Planck regime, since an almost scale- 
invariant tensor spectrum would require a certain degree 
of fine tuning, e « —qa, which however could spoil scale- 
invariance in the scalar sector— i^i^. 

The spectrum of primordial gravitational waves is not 
directly observed as after inflation it evolves to a stochas- 
tic background figw through radiation and matter eras. 
This process can be modelled in a transfer function 
T(fc>a, so tha^ 



^ ^ 1 dpg„ 
po d In k 



= T^{k)A'^, 



(92) 



where po — 3Hq/{8ttG) is the critical energy density of 
the universe {Hq is today's Hubble parameter), and /Ogw 
is the energy density of the gravitational waves. Then^ 



riT 



1 

ln/-ln/o 



In 2.29 X 10^ 



gw 



(93) 



B. Quasi-classical regime 

If one assumes that inflation takes place far enough 
from the Planck era, one can look for solutions which 
are perturbative in the small parameter ac, = 
'^nCi'2w^\ The natural choice Eq. (|64p coincides with 
classical gravity and obviously it can fit experimental 
data; so we shall concentrate on the case ac ^ 0. Rewrit- 
ing Eq. ([73]) as 

d'^Wk + cn{a ~ l)drWk + {{2a - l)k'^ 

+n^[e ~ 2 - c{a - l)]}wk = , 



it is sufficient to consider the ansatz 



„(o) 



(1) 
- acwl 



Wk — uu^ 

and solve at zero and lowest order in 
two equations. 



(94) 



(95) 

One obtains 



(0) 



dfw)^^ + [k^ + H'{e-2 

/ , — \ c 

/A 



(0) 

.(1) 



= 0, 

+ r(T) =0, 



cHdrwf + {2k' - cn')w^°^ 



(96) 
(97) 



(98) 
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the first being the usual Mukhanov equation in general 
relativity and the second having a source term. These 
equations can be solved simultaneously and exactly (even 
at higher orders in etc, using Eq. (|73p ) with standard tech- 
niques for linear differential equations with variable coef- 
ficients; however, the result for would be uninstruc- 
tively complicated and we do not show it here. Rather, 
we consider the large-scale (fc ^ Ti.) and small-scale 
{k ^ Ti.) regimes separately, and fix normalization con- 
stants by joining asymptotic solutions at horizon cross- 
ing. Equations ([M)) and ([57]) make it happen when 



nJi- 



1 
p 



(99) 



as in the classical case. In the large-scale limit the solu- 
tion is very simple and coincides with the standard con- 
stant mode 



Wk«H = Ci{l + acC2)TP , 



(100) 



where we have ignored 0{t^ p) decaying terms (we 
assume inflationary dynamics, p < —1). From 



Eq. 



w 



(0). 



at lowest order the usual relation w)^'' drW^i^"^* 



drW 



(0) 



«(327r^P[) yields at small scales 



(0) _ 



167r4i 



-ikr 



(101) 



The solution of Eq. (l97|) is given by the general solution 
of the homogeneous equation plus a particular solution 
of the inhomogeneous one. In the large k limit {Ti? terms 
neglected in Eqs. ([97|) and (|98l) ). one can show that 



w 



(1) 



46 



,(0) 



ikr 
cp — I 



IS 



The linear (in etc) part of Eq. 



(1)< 



= *(327r4i) I ^ I • (103) 



Plugging Eqs. pI7T|) and (fTU^ in Eq. pU^]) . one can fix 

the coefficients C3 = = C4. Therefore the solution at 
small scales is 



(0) 

Wk»H = W^^^ 



1 + etc 



ikr 



' cp — 1 \ a 

while the solution at large scales is Eq. (jlOOp with 



(104) 



Ci(fc) 

C2(fc) 




Csfc^P, (106) 



where = \/p{p — l)/fc and a* are defined through 
Eq. (|99p . As expected, the LQC correction term in 
Eq. (|104p decays (as t^~'^^) and becomes asymptotically 
negligible with respect to the purely classical mode. 
The tensor amplitude is 



4«2 jL2(l+p) 



where 



The tensor index reads, for small (5pi = 0(^|'j), 

n.^2(l-Hp + cpM^ -^(; + f') . 
Up to an 0(1) factor. 



<5pi oc 



(1 - cpf 



V, 



1/3 



2c 



(107) 



(108) 



(109) 



(110) 



which can be recast in terms of and M as in Eq. (|65p . 
Because the present status of the theory is not advanced 
enough to know the details of the semiclassical state, 
inhomogeneous loop quantum cosmology does not yet 
have enough predictive power to make a unique state- 
ment about the magnitude of (5pi for a given pivot scale 
/c*. However, Eqs. (|65p and (|110p illustrate the qualita- 
tive behaviour of quantum corrections in the large class 
of semiclassical states modelized by lattice refinement. 
Whether and how these states emerge in the full theory 
is a problem which goes beyond the scope of this paper. 

Actually, there is a natural scale that may fix V^^^ on 
any time slice: namely, the particle horizon. For a large 
universe this is approximated by the Hubble radius 
(we have implicitly assumed throughout the paper that 
the fiducial volume is much larger than the Hubble vol- 
ume). In Refi^, the effect of holonomy corrections were 
considered in tensor cosmological observables. In partic- 
ular, the inflationary tensor index gets an extra contribu- 
tion ^hoi oc (£piiJ,)2. Then 5p\ ^ (^pii?*)^"" < ^hoi and 
the scenario is dominated by holonomy corrections. Tak- 
ing the grand-unification scale ~ 10^** 10^^ GeV, 
(5hoi ~ lO^^'' 10^"*, so none of the corrections is observ- 
able. The same conclusion is reached even when taking 

the largest possible inverse- volume correction. By defini- 

1/3 

tion of quasi-classical regime, Vj > £pi. From Eqs. ([FT)) 
and (|58p. one has the upper bound Sp\ < O(10~^). 



IV. CONCLUSIONS 

In summary, we showed that the near-Planckian 
regime of LQC with inverse-volume corrections typically 
(but not always) leads to a strongly blue-tilted tensor 
spectrum, which is disfavoured by observations. Separate 
studies on the holonomy-corrected tensor spectrum^ are 
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not too encouraging either, as the tensor index is blue- 
tilted even in that case^i^ (during a slow-roll phase, 
ut is nearly scale invariant, therefore the inverse- volume 
contribution always dominates). Thus, the LQC superin- 
flationary phase may not be a phenomenologically viable 
scenario to explain the observational bounds on the ten- 
sor index. 

Nevertheless, there are several caveats which we must 
mention for better estimating the robustness of these re- 
sults. The first is that the near-Planckian regime occurs 
close to the bounce, a point where the scale factor is non- 
vanishing. Therefore Eq. (j67p . which was instrumental to 
the resolution of the perturbation equations, may not be 
a good approximation of the evolution of the universe 
during this first stage of expansion. A more complicated 
ansatz for the scale factor may require numerical meth- 
ods, which should be compared with the above analysis. 
A second issue is that we used perturbation theory in a 
regime which is intrinsically nonperturbative, so that the 
former may be in fact invalid and deviation from scale 
invariance be a pathology of the linear approximation. A 
treatment beyond first-order perturbations will need to 
be developed to verify this apparent tension. 

A scale-invariant or red-tilted tensor spectrum can be 
achieved in the interval ~l/qa < P ^ ^I/C^q + I), where 
we took into account the threshold for having modes ex- 
iting rather than entering the horizon. It would be in- 
teresting to check whether these values also produce a 
scale-invariant scalar spectrum. Although a very small 
tensor-to-scalar ratio could accomodate a blue-tilted ten- 
sor index, the dependence on r in Eq. (|93p is logarith- 



mic, so only a certain degree of fine tuning would al- 
low for a large blue tilt. For instance, if r ^ 10~®, 
these above bounds would still be tit < 1. In order 
to make a safe estimate of Eq. ([M]) we should plug in 
the tensor-to-scalar ratio predicted by the model, which 
we do not yet know as we still miss the scalar spectrum. 
However, anomaly cancellation may not happen in the 
scalar sector in near-Planck regime, which may actually 
be a sign that perturbation theory fails to converge. In 
fact, anomaly cancellations have been shown only when 
inverse-volume corrections are very small (a = 1 + . . . , 
where (...) ^ 1)^, that is, precisely in the large- volume, 
quasi-classical regime. Hence, for the time being the 
study of cosmological observables and consistency equa- 
tions is not negatively affected by the eventual exclusion 
of the small-volume regime from the observable window. 
This will be subject of investigation in the near future. 

Note added: Recently, tensor perturbations in the 
small- volume regime have been studied independently in 
KeiM: Where our analyses overlap, they agree. 
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